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ABSTRACT 


Model Predictive Controller (MPC) technology has been researched and 
developed to meet varied demands of need to control industrial power plants 
and petroleum refineries. This development has paved the way for the MPC 
technology too many other fields like automotive, aerospace, food processing 
industries in this paper, primary importance has been paid to the 
development of a MPC for an identified model of Multiple Input and 
Multiple Output process. In this paper, a Four Tank System has been 
considered for generation of input-output data. This data 1.e. generated input 
output data is used for the estimation of two polynomial model, name ARX 
model (Autoregressive exogenous) model and OE (Output Error) model. 
With each of model output generated, the Fit-Rates of models are compared 
to find out most efficient model. The model equations are now considered as 
plant for developing a Model Predictive Controller (MPC). Two sets of 
results are obtained after the development of MPC and tested. One is without 


noise and one is with noise. Both sets of results were a success as the output 
signals traces step input signals after some steady oscillations in real time 
with in a very short period of time which indicated a good response time. 
The MPC developed can be applied to any polynomial model with a good 
Fit-Rate, it predicts and control the process variables automatically. 
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1, INTRODUCTION 

Model identification of a single input and single output system is obtained by minimizing the output 
error of a closed loop model using constraint optimization. The optimization program is implemented by 
using FOMCON tool in Matlab [1]. A systematic procedure was derived for handling stationary random 
problems with multiple inputs and multiple outputs where only auto-spectra are known. This method can also 
be applied to nonstationary random inputs and deterministic sine and transient problems [2]. Accurate models 
were derived from the both physical and experimental data and decentralized control was demonstrated on 
the process. The main purpose was to provide connection between the abstract control theory and the real 
world. The quadruple-tank process well suited to illustrate performance limitations in multivariable control 
design caused by right half plane zeros [3]. A mathematically clean formulation of MPC emerged which 
allowed many researchers to address problems like feasibility stability, performance in a rigorous manner. 
In nonlinear area variety of issues remained which are technically complex but have potentially significant 
practical implications for stability and performance and the computational complexity necessary to achieve 
them. Simulation results of continuous stirred tank heater (CSTR) open loop and closed loop are generated 
by using Simulink implementations Matlab [4]. MPC has the ability to anticipate future events and can take 
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control accordingly. A process model is used to predict the current values of the output variables. 
The residuals the differences between the actual and the predicted outputs, serve as the feedback signal to a 
prediction block The predictions are used in two types of MPC calculations that are performed at each 
sampling instant: set point calculations and the control calculations. The coordination of the control and the 
set-point calculations is a unique feature of MPC [5] Optimization plays a significant role in designing 
process of MPC. FARX model of thermal dynamics of buildings is identified by using least square parameter 
estimation method and Compared the FARX model with the integer order ARX model [6-8]. Identified 
FARX model is validated by using IESVE software [6]. MPC is proposed to design the energy management 
system for providing reliable and optimized production of different sources of microgrid [9]. The offline 
MPC controller based on MPP method was offered for secondary loop to allow motor current to imitate the 
reference signal from the primary controller. Optimization problem in the MPC controller was altered by 
incorporating a Laguerre Model approach to decrease the number of decision variables [10]. 

The detailed above discussion is a stimulus to our present work. Modelling and controlling of 
multivariable systems is very tedious because they are highly interactive and complex [11]. Often complex 
models are beneficial for the analysis of system behavior whereas simple models are generally used for 
design of feedback controllers [12]. Design of controller system for any process to enhance the cost effective 
modelling and controlling is a difficult task [13]. The entire paper has been divided into three main parts 
namely, Data Generation [14], System Identification and Model Predictive Control. Basically, 
data generation refers to process of generating test data [15] for purpose of estimation and testing activities. 
The data needs to be generated in adherence to the set up test environment and should be in synchronization 
with developed test cases to fulfil requirement of testing model. Finally data is generated manually by 
considering the plant equations of Quadruple Tank System. System identification uses the input and output 
signals that are measured from a system to estimate the values of adjustable parameters in a given model 
structure. Model Predictive Control, advanced method of process control that has been in use in the process 
industries in chemical plants and oil refineries. They rely on dynamic model of the process, most often many 
linear empirical models obtained by system identification. 


2. RESEARCH METHOD 

In this section, data generation from quadruple tank system is generated using System Identification 
Tool Box and some brief concepts of Model Predictive Controller(MPC) are discussed. 
2.1. Physical Model 


Quadruple tank consists of two inputs namely pump voltage v1, pump voltage v2 and two outputs 
namely tank! height, tank 2 height. Schematic diagram of quadruple tank process as shown in Figure 1. 
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Figure 1. Schematic diagram of quadruple tank process [3] 
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The quadruple tank process is represented by following equations [3]: 


dhi _ al a3 yik1 
ae ~7 Vv 2ghl +7 v 2gh3 + AA DI (1) 
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“= = (2gh2 += /2gh (2) 
—_ 2gh3 oer a is (3) 
=H = 2gh4 + a ee (4) 


Where hj=level of tank 1, Ai= cross section of tank 1, aj = cross section of outlet pipe. 

The voltage applied to Pump is vj and its flow is kiv;. The parameters yi, y2 € (0,1) are determined 
from how the valves are set prior to an experiment. The flow to Tank | 1s ylkIv | and the flow to tank4 is 
(1- yl) klvl1 and similarly for Tanks 2 and 3. The acceleration of gravity is denoted by g and the measured 
level signals are k-h; and k-hz. Out of all parameters, two inputs (vl, v2) and two outputs (hl, h2) are 
considered and plotted. During the data generation white noise unmeasured disturbance signals are added to 
the two outputs of the process. These unmeasured disturbance dynamics are assumed as the stochastic 
process of following governing models [16-17] 


0.04 


C(k) = ape Nf) (5) 
Ck) = M(t) (6) 


Training data is a dataset used to fit the parameter. Validation set can be regarded as a part of the 
training set because it is used to build model. It 1s also used for parameter selection and to avoid over-fitting. 
Out of 1000 samples 750 samples are taken for training data and 250 samples are taken for validation data. 
Specifications: 

System model parameters and initial conditions of quadruple tank assumed are as follows: 

Model parameters of quadruple tank as shown in Table 1. 


Table 1. Model Parameters of Quadruple Tank [3] 


Parameters Values 
Al 28 cm? 
A2 32 cm? 
A3 28 cm? 
A4 32 cm? 
al 0.07 cm? 
a2 0.057 cm? 
a3 0.07 cm? 
a4 0.057 cm? 
k, 0.5V/cm 

g 981 cm/sec? 
kl 3.33V/cm 
k2 3.35V/cm 
yl 0.7 

y2 0.6 


2.2. System Identification 
2.2.1 Output Error Model (OE model): Structure of Output Error Model is Given by Following 
Equations [18-20]: 

y(k)=B (q-1)/F (q-1) u(k) + e(k) (7) 


1.e., non-linear optimization problem for parameter estimation of output error defined. 


¥ (V'(1) v(k) = Lies [0(K)] (8) 
Where v (k) =y(k)- x (k) (9) 
X (k) =-a1X (k-1)-a2x (k-2) +b1x (k-2) + bz x (k-3) (10) 
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2.2.2 Auto Regressive Exogenous Model (ARX model): The Equations for ARX [18, 21, 22] Model 
Given by Following: 


y(k)=B (q-1)/A (q-1) uk) + e(k) /A (q-1) (11) 
6= Mai, b1, a2, b2 Dit e’(k)) = Omin (et *e) (12) 
yw =ele = [Y-Q0]? *[Y-Q0] (13) 


2.2.3 Model Predictive Controller (MPC): The Conventional Kalman Filter [23] based on MPC 
Formulation is Developed Starting from Model: 


x(k+1) =® x(k)+P u(k)+w(k) (14) 
y(k) =C x(k)+ v(k) (15) 
& (k+j+1)k) =O% (k+j]k) +Pu(ktj]k) +Le (k+]k) (16) 

€ (k+j+1|k) = € (k+j|k) (17) 
9(k+j) =Cx (k+j]k) + € (k+jlk) (18) 


For j=0,1,2...p-1. Here, ¢ belongs to R behaves as an integrated white noise signal that provides 
compensation for the unmeasured disturbance/model plant mismatch. The initial value at beginning of 
horizon, e(k|k), 1s estimated by filtering the innovation sequence {e(k)} through a unity gain robustness filter 
[24] as follows: 


ek) = Be e(k-1) +[1-Ge] e(k) (19) 
e(k|k) =ef(k) (20) 
e(k) =y(k)-Cx(kIk-1) (21) 
& (k|k-1) = D& (k-1k-2) + Pu(k-1) + Le(k-1) (22) 


To compare the performances of MPC controllers developed using identified models, 
the performance index (JMPC) defined as: 


Juec = YS, [r(k) — y(kK)] 7 Welr(k)-y(k)] +Au(k)? Wau Au(k) [7, 25] (23) 


Is used which is based on objective function used in the proposed MPC formulation. In addition, 
the optimization technique used 1s integral square of errors (ISE) defined as: 


ISE (yi) = DNS, fyi(k) — ri(k)P [7] (24) 


Where Ns is the length of simulation trail, r(k)-y(k) is an objective function, Wg is input weight, Wau is the 
output weight. 


3. RESULTS 
3.1. Data Generation 

In this section generation of stable response for quadruple tank system identification of efficient 
model among OE model & ARX model and generation of responses for MPC with and without noise are 
obtained. Input-output data is generated from MATLAB code of four tank system without noise which is 
represented in graphical forms as shown in Figure 2a and 3b. Initial conditions of outputs such as height of 
the tank | and tank 2 given as 12.75, 13.1 respectively. 

Input-output data is generated from MATLAB code of quadruple tank system with noise which 1s 
represented in graphical forms as shown in Figure 3a and 4b. Out of 1000 samples 750 samples are taken for 
training data and 250 samples are taken for validation data as shown in above figure. 
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Figure 2. Input and outputs plots for quadruple tank system without noise 
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Figure 3. Input and output plots for quadruple tank system without noise 
3.2. System Identification 
Using data from Figure 2a and 2b, two polynomial models are generated namely OE and ARX 


model as shown in Figure 4a and 4b. 
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Measured and simulated modet ou! 




















(a) 
Figure 4. Fit rates for output yl and y2 without noise 
Model output yl has fitrates of 96.88, 71.72 for OE & ARX model and model output y2 has fit rates 


of 96.78, 73.36 for OE and ARX without noise models. Using data from Figure 4a and 4b. two polynomial 
models are generated namely ARX and OE model as shown in Figure 5a, 5b, 6a and 6b. 
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Figure 5. ARX model output yl and y2 with noise 


Discrete ARX model of four tank system 1s obtained as 


0.19081z°-14+0.1435z-2 =0.006292z7' + 0.06521z~7 


1 — 0.4197z-! — 0.409z~2 —0.05776z~—1 — 0.0187z~2 
—0.034z~1 + 0.04084z~2 0.4z°>14+ 0.11z~7 


1 — 0.4207z~-1 — 0.4367z-*. —0.0674z71 — 0.00412z~? 
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Figure 6. OE model output yl and y2 with noise 
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Discrete OE model of four tank system 1s obtained as 


0.1081z71 + 0.256z~7 —0.0817z°!1 4+ 0.144z~2 
1 —0.141z-! —0.707z-2 1—1.3580z7! + 0.3966z~2 
0.0017z~! + 0.01047z~2 0.14927°14+ 0.1397z~2 
1 — 1.78z71 + 0.7932z-2 1 — 1.88z—1 + 0.8857z-2 


From the above Figure 5, model output yl has fit rates of 86.43, 82.49 and similarly from the 
Figure 6, model output y2 has fit rates of 85.67, 84.09 for ARX and OE noise models. 


3.3. Closed Loop Response 

Simulating the MPC plant model with Simulink input and output plots are generated for MPC plant 
without noise as shown in Figure 7 and Figure 8. Prediction and control horizons are considered as 5 and 2 
respectively. Simplex optimization technique is used for optimization of MPC parameters. 
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Figure 7. Input plots for MPC plant without noise Figure 8. Output plots for MPC plant without noise 


From the Figure 8, outputs increase slowly with each step change i.e. outputs track the step inputs 
and reaches steady state. Simulating the MPC plant model with Simulink input and output plots are generated 
for MPC plant with noise as shown in Figure 9 and 10. From Figure 10, initially there is an inverse response 
for each step change and reaches steady state at the end. 
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Figure 9. Input plots for MPC plant with noise Figure 10. Output plots for MPC plant with noise 


4. CONCLUSIONS 
In this section, primary importance has been paid to the development of a MPC for an identified 


model of a Multi-Input and Multiple-Output process . A Four Tank System was considered for the generation 
of input-output data. This data, once generated, used for the estimation of two polynomial models, namely 
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ARX (Auto Regressive Exogenous) model and OE (Output-Error) model. With each of model output 
generated, the Fit-Rates of the models were compared to find out the most efficient model. After choosing 
the best model, the model equations developed. These model equations were then considered as the plant 
equations for developing a MPC Two sets of results were obtained after the MPC was developed and tested. 
One was without noise and other was with noise. Both sets of results were a success as the output signals 
traced the step input signals after some steady oscillations in real-time within a very short period of time, 
which indicated a good response time. Quite commonly MPC applications in industries involves dozens of 
inputs and outputs. The need is to develop multivariable models of such sizes through plant tests puts 
unprecedented demands on model identification techniques. In most cases, model fitting done using SISO, 
MISO methods. Because the model for each output is fitted separately in these methods, correlations among 
different outputs cannot be captured or exploited. Particularly designing a MPC in some applications where 
some of the controlled variable are either not measured or measured with large delays and must be inferred 
from secondary process measurements for satisfactory control. Many MPC's that perform well when first 
commissioned deteriorate over time, some leading to eventual shutdown. In order to sustain benefits of 
MPC's over long period of time, a mechanism to detect abnormality and dragonize its root cause is needed. 
The results can be communicated to engineers and can also be used to adapt the control parameters. 
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